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ABSTRACT 

We present our analyses of 15 months of Kepler data on KIC 10139564. We detected 57 pe- 
riodicities with a variety of properties not previously observed all together in one pulsating 
subdwarf B star. Ten of the periodicities were found in the low-frequency region, and we asso- 
ciate them with nonradial g-modes. The other periodicities were found in the high-frequency 
region, which are likely p-modes. We discovered that most of the periodicities are compo- 
nents of multiplets with a common spacing. Assuming that multiplets are caused by rotation, 
we derive a rotation period of 25.6±1.8days. The multiplets also allow us to identify the 
pulsations to an unprecedented extent for this class of pulsator. We also detect I ^ 2 mul- 
tiplets, which are sensitive to the pulsation inclination and can constrain limb darkening via 
geometric cancellation factors. While most periodicities are stable, we detected several re- 
gions that show complex patterns. Detailed analyses showed these regions are complicated 
by several factors. Two are combination frequencies that originate in the superNyquist region 
and were found to be reflected below the Nyquist frequency. The Fourier peaks are clear in 
the superNyquist region, but the orbital motion of Kepler smears the Nyquist frequency in the 
barycentric reference frame and this effect is passed on to the subNyquist reflections. Others 
are likely multiplets but unstable in amplitudes and/or frequencies. The density of periodici- 
ties also make KIC 10139564 challenging to explain using published models. This menagerie 
of properties should provide tight constraints on structural models, making this subdwarf B 
star the most promising for applying asteroseismology. 

To support our photometric analysis we have obtained spectroscopic radial-velocity mea- 
surements of KIC 10139564 using low-resolution spectra in the Balmer-line region. We did 
not find any radial-velocity variation. We used our high S/N average spectrum to improve the 
atmospheric parameters of the sdB star, deriving T e g = 3 1,859 K and log g = 5.673 dex. 

Key words: asteroseismology, stars: subdwarfs, oscillations (including pulsations); space 
vehicles: instruments. 



* Based also on observations made with the Nordic Optical Telescope, op- erated on the island of La Palma jointly by Denmark, Finland, Iceland, Nor- 
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1 INTRODUCTION 

Hot subdwarf B stars (hereafter: sdB) are horizontal branch stars 
that consist of helium-rich cores surrounded by thin hydrogen en- 
velopes. The helium core in each sdB star sustains nuclear burning 
of helium into carbon and oxygen, though in some cases the core 
may have exhausted helium at the center. Their progenitors are low- 
mass, main sequence stars, and each is likely to have undergone a 
core helium flash (M < 2Mq). Once hydrogen in the core is ex- 
hausted, it evolves towards the tip of the red giant branch, where 
an sdB progenitor retains less than 0.01 Mq of hydrogen during (or 
prior to) the helium flash. Those red giant stars lie on the high T e g 
(~30 000K) end of the horizontal branch. This small fraction of hy- 
drogen has an influence on the future evolution of sdB stars. Since 
hydrogen-shell burning cannot be sustained by a thin hydrogen en- 
velope, the stars will move directly to the white dwarf cooling track 
instead of going through the asymptotic giant branch phase. Sev- 
eral mechanisms for producing this small residual hydrogen enve- 
lope inrol\dn^_smgle-star or binary evolution have been proposed; 
e.g. iD'Cruz et all Jl996h : iHan et all d2002l 120031) . The sdB stars 
have effective temperatures and surface gravities near 30 000 K and 
log g ~ 5.5, respectively. 

It has been more than 14 year s since the discovery of the first 
pulsating sdB star (hereafter sdBV. lKilkennv et al.lll997l) . This star 
and several other stars discovered later were dominated by short 
period os cillations, which were identified as pressure (p-modes) 
modes by Charp inet et al.l dl99?t) . They are driven in the envelopes 
and have allowed us to use astero seismology to study the interiors 
of sdB stars. Several years later, iGreen et al.l (12003b reported the 
discovery of another kind of variation in sdB stars. These bright- 
ness changes were smaller i n amplitude and long er in period. Ac- 
cording to the work done bv lFontaine e t al. (2003), the new, longer 
period variations were assigned to gravity modes (g-modes) which 
dominate cooler, lower-gravity sdB stars. 

The g-modes in sdB stars probe deeper than p-modes and 
should allow us to study different parts of sdB stars. Recently, 
iReed et alj J201 l|) used g-mode period spacings for mode identi- 
fications in sdB stars. It was found that I = 1 period spacings were 
roughly 250 sec while 1 = 2 spacings were approximately 150 sec. 
There are sdB stars pulsating in both p- and g-modes at the same 
time. These stars are the most favorable, since they allow us to 
probe larger parts of those stars. At first, only p-mode dominated, 
hybrid sdBV stars were detected from the ground, but the Kepler 
spacecraft has discovered several hybrid stars w here the g-modes 
dominate across th e entire amplitude spectrum (Ree d et al1l2010l ; 
iBaranetal .11201 lij) . 

Nonradial pulsations can be described using spherical har- 
monics with quantized indices n, I, and m, where n corresponds 
to the radial order of the mode, I is the number of nodal lines on 
the surface, and m is the number of thos e lines passing through the 
poles. Following the formalism given in lLedouxl ( fl95ll) . 

0"n,;,m = O"n,i,0 + AcTu^m = 0"n,i,O + mfi(l — C n ,l) (1) 

where <7 n ,;,m is the shifted frequency and a n ,i,o is the (degener- 
ate) frequency without rotation. If the Ledoux constant C ni ; is es- 
timated and the multiplet splittings Aa Ut i >m are measured, we can 
derive the rotation frequency £1. It should be remembered that this 
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Table 1. Details on KIC 10139564. The last row reports the number of 
dates used in our analysis. Starting and ending points are given in BJD and 
are shifted to 2455000.0 epoch. 



R.A. [2000.0] 


DEC. [2000.0] 


19 24 58.16 


+47 07 53.6 




log g [dex] 


log N He /N H 


16.1 31,859(126) 


5.673(26) 


-2.201(36) 


start end 


length [days] 


points 


276.485140 738.935300 


~462.5 


625974 



equation is valid for slow, solid-body rotators with pulsation pe- 
riods much shorter than the rotation period. In the absence of ro- 
tation, all modes of the same degree but with differing m values 
are present in pulsating stars, but they are degenerate in frequency, 
though they differ in angular eigenfunctions. In this case, we only 
see the a n .ifi frequency. However, rotation lifts the degeneracy in 
m, resulting in (21 + 1) periodicities with different frequencies. 
Periodicities with positive and negative m are shifted in the oppo- 
site direction of <Jn,i,o- The m = frequency remains unaffected by 
rotation. 

It is generally accepted that p-modes detected in sdB stars are 
of low order and of low degree. The Ledoux constant for such 
modes is typically close to zero. However, for the fundamental 
and first overtone I = 2 , 3 modes, C n> i may have values up to 0.3 
dCharpinet et~aUl2002l) . This increase in the Ledoux constant cre- 
ates a decrease in the frequency splittings. In the case of high-order 
g-modes, the Ledoux constant depends on the I parameter accord- 
ing to the following equation: C n ,; < (I 2 + ■ It is around 0.5 
for I = 1 and close to 0.16 for 1 = 2. If we assume a star rotates as 
a solid body, the ratio of the Ledoux constants between the p- and 
g-mode regions can be used to constrain the degree of the g-modes. 
A p/g splitting ratio of 0.55 indicates 1=1 g-modes and a ratio of 
0.81 indicates I = 2 g-modes. 

Observations of long period sdBV stars, including detection 
of new pulsators, are not easy from the ground. Relatively long 
period photometric variations (with only a few cycles per night) 
can be affected by variable sky transparency. To make a defini- 
tive determination of the pulsation frequencies requires extended 
photometric campaigns, preferably at several sites widely spaced 
in longitude to reduce day/night aliasing. However, Kepler pro- 
vides long-duration, continuous, homogenous, evenly spaced time- 
series photometry, making it an ideal instrument for asteroseismol- 
ogy. The Kepler science goals, mission design, and overall per- 
for mance are reviewed b y iBorucki et al.l d2010l) , iKoch et all d20ld) 
and I Jenkins et al.l d2010h . This is currently the state-of-the-art way 
to observe variable stars. With this instrument, we can reach down 
to 19 mag and still achieve better signal-to-noise (S/N) observations 
than from the ground. 

In this paper, we present analyses on KIC 10139564 
(2MASS 19245816+4707536), a p-mode dominated sdBV star de- 
tected in Kepler data. Our analyses are completed on combined Q5 
through Q9 data spanning 15 months of continuous observations. 
Such long, continuous coverage - never achieved before for com- 
pact pulsators - significantly improves the S/N ratio in the ampli- 
tude spectra. This coverage allowed the detection of more frequen- 
cies and helped to improve mode identification, as c ompared to the 
prelim inary results from the one-month survey data (Kawa ler et al.l 
|2010|) . In addition to Kepler data, we collected spectroscopic data, 
which we present in the next section. 
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2 SPECTROSCOPIC DATA 

On the nights of 31 May (2 spectra), 19 June (2) and 23 July 2010 
(6), we obtained 10 spectra in total of KIC 10139564 to monitor 
possible radial-velocity variations and, thus, check its binary status. 
We used ALFOSC at the Nordic Optical Telescope with grism #16 
that yield spectra of the H^-He region with dispersion of 0.77A/pix 
and a resolution of 2.2A, with a median S/N = 48 per pixel. The 
exposure times were 900 seconds. The spectra of 23 July were not 
taken consecutively so cover 0.3 days. 

The data were homogeneously reduced and analyzed. Stan- 
dard reduction steps within IRAF were bias subtraction, removal of 
pixel-to-pixel sensitivity variations, optimal spectral extraction, and 
wavelength calibration based on arc lamp spectra. The target spec- 
tra and mid-exposure times were shifted to the barycentric frame of 
the solar system. The spectra were normalized to place the contin- 
uum at unity by comparing with a model spectrum for a star with 
similar physical parameters as we find for KIC 10139564. 

Radial velocities were derived with the FXCOR package in 
IRAF, with resulting errors of about 10 km/s per spectrum. We used 
the lines of Hp, H 7 , and to determine the radial velocities. 

We find an average radial velocity of 1.0 ±2.9 km/s and find 
all data points except the two from 31 May consistent within la 
with average radial velocity. The two measurements of 31 May 
combine to a radial velocity of 15.5 ±7.1 km/s. We conclude that 
there is no clear evidence for orbital radial- velocity variations for 
periods shorter than half a day. 

The 10 individual spectra were combined into an average 
spectrum with S/N~150, which was used for an improved deter- 
mination of the atmospheric parameters of KIC 10139564. We fit- 
ted the mean s pectrum using the same method and grid of models 
as that used in 0stensen et al. (2010a). (For a description of these 
models, see iHeber. Reid &*W erner 2000). The fit is shown in Fig- 
ureQ] The physical parameters resulting from the fitting procedure 
are listed in TableQ] and places KIC 10139564 at a temperature 
some 600 K lower and at a surfac e gravity 0. 14 dex lower than that 
obtained by 0stensen et alj d2010al) whose data were based on a 
single low S/N spectrum. 



-20 -10 10 20 
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Figure 1. The fit of an LTE model spectrum to a mean of 10 spectra obtained 
with the NOT/ ALFOSC. Note that the H e line was kept out of the fit due to 
interstellar Ca II contamination. 

3 PHOTOMETRIC DATA 

All data included in this paper were obtained through the Kepler 
Asteroseismic Science Consortium (KASCfl Since pixel data for 
Q6 and Q8 are still not available to us, we used flux measurements 
already calculated using default apertu res defined by the in-house 
extraction pipeline jjenkins et al.12 010). We have chosen to use the 
raw fluxes, which are calibrated for instrumental effects (bias, flat- 
field, cosmic rays) but not corrected for on-board systematics or 
flux excesses (often called contamination). Contamination does not 
influence our main goal, since we rely on periods and not ampli- 
tudes. We have corrected the data for on-board systematics to re- 
move discontinuities in the data. 

TableQ] contains an observational report and the number of 
points used in our analysis. It is worth noting the unprecedented 
length of continuous coverage which cannot be obtained from the 
ground. The duty cycle is 92%. We decided not to plot the pho- 
tometric data, sinc e they do not differ f rom those presented in the 
previous paper by iKawaler eTall feOlCh other than being substan- 
tially longer in time. 

We analyzed short cadence (SC) data. Each data point is a sum 
of nine integrations of 6.02 s followed by a 0.52 s readout, resulting 
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in an SC sampling of 58.85 s. This translates into a Nyquist fre- 
quency of 8496.356 /iHz in the spacecraft reference frame. Time 
stamps given in that frame were corrected to the barycenter of the 
Solar System. This correction has smeared the Nyquist frequency 
because data are no longer evenly sampled. 

In preparing data for Fourier analysis, we first removed all 
SC points flagged with errant times or fluxes. Then, all remaining 
points were subject to detrending which was done via a cubic spline 
fit calculated on 0.5 day subsets and subtracted from the original 
data. Detrending with this interval will strongly suppress all fre- 
quencies below 12 /jHz. Then we removed all outliers using 3.5<r 
clipping before the data were subjected to Fourier analysis. We did 
not include the 30 days of data from the survey phase, since they 
were too temporally separated from the data presented here, which 
could cause a complex pattern of peaks in the amplitude spectra 
and hinder our prewhitening process. 

During Q5, Q6, Q7 and Q9, there were no safe-mode events, 
and cadences were only lost while the spacecraft was not in fine 
pointing (which occurred during monthly data transmissions to 
Earth). This makes these four quarters failure free. Unfortunately 
during Q8, the spacecraft experienced three safe-mode events: one 
which delayed the start of Q8 by about two weeks, one which 
caused a three-day gap, and one which prematurely ended Q8 ob- 
servations 6 days early. 



4 AMPLITUDE SPECTRUM 

We used Fourier analysis to identify candidate pulsation frequen- 
cies. Then we completed a nonlinear least-squares fit, including 
each periodicity in the form Ai sin(ujit + (f>i). This iterative process 
followed the standard prewhitening procedure and continued until 
all peaks with amplitudes above a threshold of S/N = 4 had been 
removed. While a threshold of S/N=4 is the commonly-used limit 
for sdB pulsations, it is not a strict limit and so frequencies near 
this limit (S/N^4.5) may be considered tentative, and are marked 
as such in Table|2] Continued observations by Kepler may be used 
to confirm or refute them. The frequency resolution is 0.0375 /iHz, 
defined as 1.5/T, where T is the time baseline of the data. We have 
coverage fifteen times longer than the 30 days of survey data, and 
the noise level in the amplitude spectra are lower roughly by a fac- 
tor of \/l5 = 3.87. 

When the spacecraft is rolled every 3 months, the stars move 
onto different CCDs. Moreover, any deviation in pointing accu- 
racy and other on-board systematics may cause different flux levels 
between quarters, which in turn leads to a different amplitude of 
light variation. Although we detrended the data using 0.5 day sub- 
sets and 4a clipping, it is possible that instrumental systematics 
remain in the data. Such instrumental flux variations would de- 
crease the effectiveness of prewhitening and manifest themselves 
as increased residuals and residual peaks with fork-like structures 
(multiple, closely-spaced peaks). Since many of the residuals in 
KIC 10139564 do show fork-like structures we examined other 
options to further correct the data for instrumental flux changes. 
The best possible method would be to use a source external to 
the spacecraft, such as the star itself, for systematic flux correc- 
tions. Pulsation amplitudes can be used for this purpose, so long 
as the amplitudes themselves are stable with time. While many 
sdB pulsation amplitudes a re variable, some do appear to be re- 
markably stable (see §4.2 o f lReed et alj20l3) . As such, we sought 
to correct these data using the amplitudes of the pulsations them- 
selves. While this re-normalization could affect pulsation ampli- 



tudes, we are not interested in absolute amplitudes but rather the 
best possible frequency solution, one measure of which are resid- 
uals after prewhitening. The procedure we used was to split the 
data into month-long segments and fit the highest-amplitude pul- 
sations. We then used those amplitudes to correct the original data 
on a month-by-month basis and re-prewhitened and analyzed the 
resultant residuals. We did this normalization one high amplitude 
pulsation frequency at a time. Only for the peak fig (Table[2} did 
we achieve significant improvement. An example region is shown 
in Figuref2] In this figure, the top panels show the unnormalized 
data and residuals while the bottom panels show the same for the 
amplitude-normalized data. The frequency solution is the same in 
both cases but the normalized residuals are smaller with less fork- 
like structure. The detection thresholds (horizontal dotted lines) for 
un- and normalized data are equal to 0.034 and 0.033 ppt, respec- 
tively. Since the residuals were reduced for all frequencies, any am- 
plitude changes in fig were most likely systematic and not intrinsic 
to that peak. For KIC10139564 this procedure did not change the 
frequency solution and so the sole benefit was in reducing several 
of the residual peaks below the detection threshold and removing 
some of the fork-like structures, which is a purely cosmetic result. 

For this amplitude-normalized dataset we calculated the de- 
tection threshold as four times the average level of the residual 
spectrum in the range of 230 /iHz to the Nyquist frequency. The 
detection threshold is 0.033 ppt (parts per thousand), but can devi- 
ate locally by 3% when using 500 /iHz bins in signal-free regions. 
Please note that observed amplitudes are smeared by a finite sam- 
pling, leading to a decrease of 15-20% for most periods and up to 
35% for the shortest periods. 

KIC 10139564 is the only short period dominated (p-mode) 
pulsator among sdBV stars detected by Kepler. Based on 30 days 
of survey data. lKawaler et al . (2010) noted that the amplitude spec- 
trum of this star is dominated by several peaks between 160 and 
200 s (6200 and 5050 /iHz) with the highest amplitude being about 
9.9 ppt. One single peak in the low-frequency region was also men- 
tioned, but it had yet to be confirmed. Our analysis confirms all 
signals reported in the previous analysis. However, our frequency 
resolution is fifteen times better than in the survey data, allowing 
us to better resolve crowded peaks that lead to small changes in fre- 
quencies. In total, we have identified and successfully fit 57 period- 
icities with amplitudes above our detection threshold of 0.033 ppt 
(Tabled ■ Note that some signals prewhitened in the previous anal- 
ysis were not prewhitened in ours. With our improved resolution, 
we found that the signal in some regions is very complex and diffi- 
cult to fit. We discuss these regions in E|5 . 5 1 
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Table 2. The frequency list as a result of our prewhitening process for KIC 10139564. The detection threshold is 0.033 ppt. The numbers in parentheses are the 
1-cr fitting errors of the last digits. This list contains only peaks prewhitened from our data, which does not include some complex signal patterns described in 
£|5.5| and LC artifacts. Observed amplitudes are smeared by a finite sampling, leading to a decrease of 15-20% in most periods and 35% for the shortest periods. 
The asterisks denote tentative peaks with S/N^4.5. The single horizontal line separates low and high-frequency regions. Columns six and seven provide I and 
m identification as a result of our analys is. This identification is based on full or partial multiplets. The last column, denoted K10, provides the notation of 
frequencies used in iKawaleretallfeOlOl) . 



ID 


Frequency [/iHz] 


Period[s] 


Amplitude [ppt] 


S/N 


/ 


m 


K10 


f A 


315.5766(21) 


3168.803(21) 


0.050(8) 


6.0 




-1 




f B 


315.8224(8) 


3166.336(7) 


0.138(8) 


16.6 


1 







fc 


316.06559(31) 


3163.9002(31) 


0.339(8) 


40.8 




+1 


fO 


f D 


345.2455(28) 


2896.489(24) 


0.038(8) 


4.5 


2 







I, 


345.9747(31) 


2890.385(26) 


0.034(8) 


4.1 




+2 




If 


394.0282(14) 


2537.889(9) 


0.078(8) 


9.4 


1 







f G 


394.2876(14) 


2536.220(9) 


0.079(8) 


9.5 




+1 




% 


518.8984(19) 


1927.160(7) 


0.057(8) 


6.8 




-1 




fi 


519.1540(14) 


1926.2106(50) 


0.079(8) 


9.4 


1 







fj 


519.4022(13) 


1925.2904(49) 


0.081(8) 


9.7 




+1 




f? 


2212.6076(29) 


451.95543(6) 


0.037(8) 


4.4 










3276.9340(30) 


305.163302(28) 


0.035(8) 


4.2 








h 


3540.4565(22) 


282.44945(18) 


0.048(8) 


5.8 


2 







h 


3541.4320(19) 


282.37165(15) 


0.056(8) 


6.8 




+2 




h 


5048.7427(21) 


198.06912(8) 


0.052(8) 


6.2 




-2 




h 


5049.7109(23) 


198.03114(9) 


0.046(8) 


5.5 


2 







ti 


5052.6025(31) 


197.91781(12) 


0.034(8) 


4.1 








fa 


5286.14933(22) 


189.173619(8) 


0.478(8) 


57.4 




-2 


fl 


fg 


5286.56176(20) 


189.158861(7) 


0.544(8) 


65.4 




-1 




fio 


5286.97660(15) 


189.144018(5) 


0.723(8) 


86.9 


2 







fn 


5287.39188(7) 


189.1291628(25) 


1.521(8) 


182.9 




+ 1 


f2 


fl2 


5287.80529(12) 


189.1143764(42) 


0.899(8) 


108.1 




+2 


f3 


fl3 


5410.6980(9) 


184.819039(30) 


0.121(8) 


14.6 




-2 




fl4 


5411.60111(50) 


184.788195(17) 


0.210(8) 


25.3 


2 







fl5 


5412.5138(7) 


184.757036(25) 


0.143(8) 


17.2 




+2 




fie 


5413.39487(29) 


184.726964(10) 


0.363(8) 


43.6 


1 





f4 


fir 


5413.8205(9) 


184.712442(29) 


0.124(8) 


14.9 




+1 




fl8 


5472.86225(3) 


182.7197459(10) 


3.712(8) 


446.4 








f7 


fl9 


5570.0278(10) 


179.532318(31) 


0.111(8) 


13.4 




-3 






5570.4811(31) 


179.51771(10) 


0.034(8) 


4.1 




-2 




f 2 l 


5571.3956(15) 


179.488243(49) 


0.070(8) 


8.4 


3 







f22 


5572.2929(23) 


179.45934(8) 


0.045(8) 


5.5 




+2 






5572.7357(29) 


179.44508(9) 


0.037(8) 


4.5 




+3 




f'24 


5616.4852(26) 


178.04730(8) 


0.040(8) 


4.8 








Q 5 


5618.9672(29) 


177.96865(9) 


0.036(8) 


4.3 








t'2(i 


5708.5701(21) 


175.17522(6) 


0.051(8) 


6.2 








f27 


5709.3992(28) 


175.14978(8) 


0.039(8) 


4.6 








f28 


5746.6089(13) 


174.015672(41) 


0.079(8) 


9.5 




-3 




f29 


5746.9531(26) 


174.00525(8) 


0.041(8) 


4.9 




- 




fill 


5747.0936(6) 


174.000993(17) 


0.186(8) 


22.3 




-2 




f 3 l 


5748.0723(8) 


173.971367(23) 


0.142(8) 


17.1 


3 







f32 


5749.0647(12) 


173.941338(35) 


0.091(8) 


10.9 




+2 




f33 


5760.16218(1) 


173.6062229(4) 


7.950(8) 


955.8 




-1 


f8 


f34 


5760.58982(4) 


173.5933352(12) 


2.741(8) 
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Figure 2. Top panels show the unnormalized data and residuals while the bottom panels show the same for the amplitude-normalized data. The dashed 
horizontal line (here and in all other figures) represents the detection threshold of 0.033 ppt for normalized and 0.034 ppt for unnormalized data. 
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Figure 3. Top panel: amplitude spectrum in the low-frequency region (note 
that abscissa axes are expressed in period instead of frequency and the top 
panel is truncated in amplitude). Three groups of peaks are evident. Num- 
bers in the plots (here and in all other figures) denote separations between 
peaks, and they are given in abscissa axis units. Bottom panel: residual spec- 
trum after removal of nine peaks (f a to f j) in this region. 



dotted line, between them could create a multiplet with a splitting 
of 0.36 fiHz on average. This value is expected for a quintuplet 
since the Ledoux constant is 67% smaller for I = 2 modes. These 
frequencies also overlap an artifact listed in Data Release Notes 12 
Notes as "wide" and it is possible they are not intrinsic to the star. 
However, that artifact was discovered by one of us (ASB) in Q0 
data. In addition, our two peaks do not resemble the artifact found 
in Q0, so we infer them to be produced by the star. 

All peaks except fo and fE are organized into three groups 
with the separations between them of 610 and 630 sec, respectively. 
This is not suiprising, since g-modes should be equally spaced in 
period in accordance with asymptotic theory. In KIC 10139564 we 
could have a sequence of I = 2 modes but only with every fourth 
multiple (620/4=155) detected. However we assigned them as I = 1 
modes based on their multiplet structures, and that makes them a 
sequence with spacings of (620/2=)310 or (620/3=)207 sec, if only 
the second or third multiples were detected, respectively. While the 
I = 1 period spacings would not be consistent with other pulsating 
sdB stars, those stars are g-mode dominated while KIC 10139564 
is p-mode dominated and the frequency splittings indicate these g- 
modes are likely 1 = 1. 



5 DISCUSSION 

5.1 Signal in the low-frequency region 

In the low-frequency region, we have detected ten peaks (fA to fj). 
All low-frequency peaks are listed in the upper part of Tableland 
are shown in Figure|3] We associated all of them with g-modes. 

In each region, the signal is split into two or three significant 
peaks as seen in Figure|4] The frequency separations between peaks 
are similar and includes periodicities with amplitudes that have 
S/N<4, shown in Figure[4] In Figure|5] we plot the window func- 
tion to show that the splitting is not a result of the aliasing pattern. 
Since rotation can remove degeneracies in pulsation frequencies of 
modes with the same I but different m value, we suggest that these 
are rotationally split multiplets. With two or three peaks observed, 
we can consider them to be I = 1 triplets; however with more data, 
it is possible they may turn into quintuplets. Two frequencies fu 
and fE are different. They show larger splittings of about 0.72 pHz. 
An insignificant frequency (S/N=3.8), marked with a vertical dash- 



5.2 Signal in the high-frequency region 

The high-frequency region contains 47 significant peaks. We show 
this region in Figure[6] The highest peak has an amplitude of 
7.950 ppt and four other peaks are between one and six ppt. Most 
of the peaks, however, are of low amplitudes, barely exceeding 
S/N= 10. Signals prewhitened from our data are denoted as p, LC 
means long cadence artifacts and s stands for complex signal re- 
gions, which we were not able to remove from the data using 
standard techniques. LC artifacts are peaks at the long cadence 
frequency (566.41 /iHz) and its harmonics up to the Nyquist fre- 
quency. They were identified during the first year of the Kepler 
mission, and their frequencies are given in the Kepler Instrument 
Handbook (KIH). Besides LC artifacts, Kepler data may also con- 
tain other spurious peaks. Some of them have been identified in 
KIH, but there are others that are still not documented very well. 
iBaran et al] d201ll) reported another set of artifact signals detected 
in many stars analyzing Kepler public data. It is also possible that 
artifacts may come and go, and a frequency of a specific artifact 
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Figure 4. Zoomed amplitude spectra for the g-mode region. Each pair 
of panels contains the original (top) and residual spectra (bottom). 
Prewhitenned peaks are marked with arrows. Frequency separations, which 
include two members with S/N<4 (marked with vertical dash-dotted lines), 
are labeled in the residual panels; those separations are not precise and may 
differ by up to 0.03 /iHz. 



Figure 5. The window function calculated from Q5-Q9 data. 
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Figure 6. The high-frequency region of KIC 10139564. This plot is trun- 
cated at 0.5 ppt to show low-amplitude peaks. Labels help to see how this 
region is organized. We found at least one strong LC artifact, many peaks 
we could prewhiten (p) and complex signal (s) which we left untouched. 



may also change. In the amplitude spectra we have detected all LC 
harmonics starting from the 6 th up to the 14 th , which is the Nyquist 
frequency. 

The first few periodicities (fi to f4) are well separated from 
the bulk of p-modes. Frequencies fi and f2 appeared to be singlets. 
At 3540 /iHz, we fit two significant peaks (f3 and f4) as seen in 
Figure|7] In addition, two other peaks seem likely to represent real 
signal (particularly the peak close to 3541 /iHz with S/N = 3.6), but 
we need more data to confirm this. Both significant peaks have low 
amplitudes and the separation between them is 0.9755 /iHz. Esti- 
mating frequencies for the two other peaks and calculating sepa- 
rations reveals that all four peaks create an almost equally spaced 
quadruplet. We found a similar case at 5049 /iHz shown in Figure[8] 
We detected two low-amplitude, yet significant, peaks (fs and fe) 
separated by 0.9682 /iHz with another insignificant (S/N=3) peak 
0.482 /iHz away from fe. The separations between peaks in the 
group at 3540 and 5049 /iHz are similar. The peak at f7 stands alone 
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Figure 7. The frequency region around 3540 /iHz. We detected two sig- 
nificant peaks f3 and f4, but there are two other peaks slightly below the 
detection limit which fit into the evenly spaced structure. 
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Figure 8. The frequency region around 5049 fiHz. Similar to Figure[7] we 
found two low-amplitude peaks fs and f6 and another peak slightly be- 
low the significance limit. It is worth noting that separations between these 
peaks are comparable to the peaks in the Figure[7] 
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Figure 10. The frequency region around 5412 fiHz. This region shows five 
peaks (fi3 to fi7), but unlike the previous region, they do not create an 
evenly spaced obvious quintuplet. 
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Figure 9. The frequency region around 5287 fiHz. This region contains a 
clearly detected quintuplet (fg to fi2). 



and barely exceeds the detection threshold. If f/7 is a part of a mul- 
tiplet, other components are currently undetected. 

The region presented in Figure|9] contains a clearly resolved 
quintuplet (fg to fi2) with separations of ~0.41 fiHz, comparable 
to those in the previous p-mode regions. As may be seen in the 
residual spectrum (bottom panel), all peaks prewhitened with al- 
most no residual signal left. 

Figure[l0] shows the region near 5412 fiHz that also contains 
five significant peaks (fi3 to fi7). However, these peaks do not cre- 
ate an evenly spaced structure, yet the separations look familiar. 
One separation is ~ 0.4256 /iHz, while the other three appear ap- 
proximately twice that near 0.9 /jHz. In addition, after all signifi- 
cant peaks were removed, we found another peak (with S/N=3.9) 
between fi3 and fi4 that may fit into this sequence. Two possible 
interpretations exist. Either all the peaks create one large / >3 mul- 
tiplet with separations of ~0.44 /iHz or, as we prefer, fi3 to /ig are 
part of a quintuplet, fi6 and fi7 are part of a triplet and the two 
multiplets just happen to be separated by ~0.9 /iHz. 

The peak fi8 shown in Figure[TT]is different. It seems to be 
single with some low-amplitude complex signal at slightly lower 
frequencies. The amplitude of this peak was used to normalize data 
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Figure 11. The frequency region around 5473 /jHz with only one peak fig 
and a wide, low amplitude signal at lower frequency. 
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Figure 12. The frequency region around 5572 fiHz. We detected five peaks 
fl9 to f23 with separations comparable to the previous regions. This may 
represent a septuplet. 
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Figure 13. The frequency region near 5748 /tHz with five peaks (f28 to f%2 )■ 
Peaks f3o to {32 might be part of a quintuplet, since the separation between 
them is twice the smallest splitting in other multiplets. 
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Figure 14. The region around 5760 /iHz which contains a triplet (f33 to 
f35) with splittings comparable to other regions. This would indicate that 
the highest amplitude peak is not a radial mode. 



between individual months, which significantly removed residuals. 
This is the only single high amplitude peak we have detected in our 
data, and its lack of multiplets, despite its high amplitude, leads us 
to conclude it is a radial mode. In the residual spectrum (bottom 
panel of FigurefTTt. we notice an obvious and very wide signal ex- 
cess. We did not try to prewhiten any peaks from this "forest". It 
resembles other regions of peaks between 5900 and 6400 /iHz. 

The signal at 5572 /iHz, shown in Figure[l2] might be of spe- 
cial interest. We detect five peaks, fig to f23, with an average spac- 
ing of 0.45 /iHz assuming that two peaks just below the detection 
threshold are intrinsic to KIC 10139564. The residual spectrum has 
a peak at the expected frequency to complete the sequence but it is 
not significant (S/N=3.41). This region could contain an 1 = 3 sep- 
tupled Thus far, the detection of 1>2 modes was generally regarded 
as unlikely, since these modes suffer from large surface cancellation 
effects, making their observed amplitudes greatly reduced, assum- 
ing the intrinsic amplitudes to be similar. With Kepler data we can 
reach down to the ppm level, making the detection of septuplets 
possible. More data are definitely needed to confirm the nature of 
the peaks in this region, but it is a good candidate for the first de- 
tection of I ^ 3 modes in any sdB star. 
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Figure 15. The amplitude spectra of the highest amplitude periodicities 0^33 
to f35) calculated from data split into quarters. Notice the amplitude shift 
between the outer frequencies. In Q8 the peak on the right is slightly higher 
than the one on the left. 



We do not show the next region at 5618 /iHz, since it only con- 
tains two low amplitude peaks f?4 and f25 separated by 2.482 /iHz. 
It may be the case of another septuplet with the two peaks being 
the first and last in a sequence of seven modes separated approx- 
imately by 0.42 /iHz. The residual spectrum shows one peak near 
the detection threshold and fits the potential sequence well. Sim- 
ilarly, only two low-amplitude peaks f26 and f27 were found at 
5709 /tHz. More data are necessary to determine the nature of these 
two regions. 

Similar to the region shown in Figure[T2] we have detected 
a triplet at 5748 /iHz (f3o to f32) with separations of 0.98 /tHz on 
average (FigurefTTt. This could easily fit a quintuplet with intrin- 
sic splittings of 0.49 /iHz, similar to other multiplets. In addition, 
f 28 and the peak (S/N~4) at 5747.55 /iHz could potentially fit this 
sequence. Another option could be that f28 to f29 are part of a dif- 
ferent multiplet with shorter separations of 0.3442 /iHz which is 
similar to the multiplets found at 5840 and 6235 /iHz. 

The next p-mode region is shown in Figure[l4] It is centered 
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Figure 16. The frequency region around 5840 /jHz. We detected two peaks 
f'36 andf37. 

around 5760 /iHz and contains three significant peaks (f33 to f3s). 
Since radial modes suffer no geometric cancellation, it is usually 
anticipated that the highest amplitude modes are, therefore, radial. 
This does not seem to be the case for KIC 10139564. All three 
peaks are spaced by the distance of roughly 0.42 /iHz, forming a 
triplet. These spacings in the highest amplitude frequencies provide 
assurance that this frequency spacing, which is common to many 
regions of KIC 10139564, is intrinsic and related to rotationally 
split multiplets. The residual spectrum shown in the bottom panel 
of Fi gure[T4l indicate s that all those peaks are unstable in their am- 
plitudes/phases. In Figure[T5"lwe show how the peaks change their 
amplitudes between quarters; in Q8 the amplitude of the right peak 
was the highest. 

The region at 5840 /iHz contains low amplitude peaks fz§ and 
f37 with very complex shapes shown in Figure[l6] They appear to 
be suffering from amplitude/phase variations or as unresolved fre- 
quencies, which makes our fitting process difficult and our estima- 
tion of their frequencies and amplitudes less precise. Regardless, 
we have fitted f36 and f37 with a separation of 2.1953 /iHz. Other 
peaks in this region cannot be fitted. We will discuss this region in 
detail in ^531 

The region at 6001 /iHz is one with complex signal patterns 
and is described in Sect |5.5l Two frequencies and f3g look rel- 
atively high and stable in amplitudes, so we fitted and prewhitened 
them (listed in Table|2j. The next two regions around 6076 /iHz 
(fzio) and 6106 /iHz (£41) are not shown since we fitted just one sig- 
nificant low amplitude peak each. The region at 6235 /iHz shows 
three significant peaks (FigurefTTt. but their structure is very com- 
plex. We could easily prewhiten two of them (£42 and £43), leaving 
the third one untouched. The residual spectrum clearly shows am- 
plitude and/or phase variation of all peaks. The smallest (lowest 
frequency) signal has a fork-like appearance, which may mean it is 
unresolved in the original spectrum, and we decided not to make 
an attempt to remove it. The spacing between the two prewhitened 
peaks is slightly smaller than the common value of 0.45 /iHz, as 
derived in other regions. This deviation can be explained by a dif- 
ference in the Ledoux constants (SjT). The remaining signal appears 
at exactly the right frequency to fit a triplet, leading us to conclude 
this region is another rotational multiplet. There is another signifi- 
cant signal in the amplitude spectrum around 7633 /iHz with a sin- 
gle frequency £44 detected. 

The region at 81 18 /iHz is the highest frequency region with 
three peaks f45 to £47 that we could easily prewhiten. In fact, it is 
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Figure 17. The frequency region around 6235 /iHz with two peaks detected 
(142 and 143). The wide peaks possibly indicate amplitude/phase variation. 
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Figure 18. The frequency region around 8118 /iHz with three peaks de- 
tected (£45 to £47). The rightmost peak is very wide and is likely unresolved, 
as the residual seems to confirm. Including the two near the detection limit, 
this region form a quintuplet. 

very close to the Nyquist frequency, leaving some doubts whether 
this region is mirrored across that frequency. If so, the real signal 
could appear at approximately 8874 /iHz instead. Similar to the first 
region at 3540 /iHz, this region is also well separated from the bulk 
of signals we detected in the amplitude spectrum. We show this 
region in Figure[T8] It seems obvious that five peaks are spaced by 
the common splitting (two just below the detection threshold), and 
we interpret them as another I = 2 quintuplet. 

5.3 Multiplets 

We expect all nonradial modes of a slowly rotating star to be split 
according to the Ledoux equation (Eq.[T}. This is indeed what we 
see in the amplitude spectrum of KIC 10139564. Most peaks are 
split with roughly consistent splitting among modes of the same 
kind (p- or g-modes). 

In Figure[T9l we show spacings between p-modes with 
S/N>5.5. Most group at approximately 0.45 /iHz, with a few hav- 
ing double that spacing at 0.9 /iHz and two with three times that 
spacing near 1.4 /iHz. We interpret the group with the smallest 
spacings to arise from modes of consecutive m values (Am = 1), 
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Figure 19. A histogram of frequency spacings in the p-mode region. 

the next group due to Am = 2, and the final pair due to Am = 3. 
Using the average splitting between p-modes of 0.452±0.030 /jHz) 
and C n .i = (splittings below 0.4 /jHz were not taken into account 
since we suspect they are caused by a non-zero Ledoux constant), 
we derived an average rotation period of 25.6±1.8 days. 

Multiplets at low frequencies have average splittings of 
0.2504±0.0019 and 0.3646±0.0042 ^Hz. As mentioned in Jj]the 
ratio of the Ledoux constant between the p- and g-mode regions 
can be used to constrain the degree of the g-modes. The smaller 
splitting has a ratio of 0.55 and the larger splitting has a ratio of 
0.81 which matches the Ledoux ratio if the g-modes are I = 1 and 
I = 2 , respectively . 

iBaran et alj (2009) reported the discovery of a symmetrical but 
not evenly spaced quintuplet in the sdBV star Balloon 090100001 
(Bal09) and explained this uneven spacing as differential rotation 
on the surface. In KIC 10139564, we detected the quintuplet fs to 
fi2, which is symmetric, but components with \m\ = 1 are spread 
wider (0.41505±0.00042 /jHz on average) compared to the |m| =2 
modes (0.41290±0.00044 ^Hz). This leads us to the conclusion 
that rotation in KIC 10139564 may be slow est on the equator which 
would be opposite to the result obtained bv lBaran et alj d2009l) . 

Pulsation amplitudes of identified multiplet components can 
also be use d to constrain the puls ation inclination (supplementary 
materials o f lCharpinet etalfeOl J). If the components of individual 
multiplets are intrinsically driven to the same amplitude, then ob- 
served amplitude differences are caused by the inclination at which 
they are viewed. This is particularly true for I 3 modes which 
have a higher degree of geometric cancellati on (Figure 8.5 of sup- 
plementary material of ICharpinet et alj|201lb . For KIC 10139564 
we detected four possible triplets (some with components below 
the detection threshold) yet they do not have consistent amplitudes 
between the various components. This is particularly evident for 
the complete quintuplet where the m = +1 and +2 components 
have triple and double the amplitudes of their counterparts, respec- 
tively. Likewise, we observe a triplet with spacings which indicate 
we are seeing the m = —2, +1 and +2 components, which can- 
not be attributed to viewing inclination. As such, the amplitudes are 
not sufficiently consistent to be used for constraining the pulsation 
inclination. 

5.4 Companions to KIC 10139564 

ISilvotti et al.1 d2007l) have detected the first planet around an sdB 
star. They showed that pulsations detected in sdBV stars can be 
used to search for periodic variations of pulsation phases, which 
could be caused by orbiting bodies. However, the pulsation phases 
have to be monitored for a minimum of 2 to 3 years and measured 
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with low uncertainties. It means that resolved high-frequency and 
high-amplitude modes are preferred. With several high-amplitude 
modes, KIC 10139564 is one of the best candidates where we can 
look for phase variations in terms of planet hunting. KIC 10139564 
has been observed for almost two years, but at the time of this anal- 
ysis, we had access to only 15 months of data. There are a few high- 
amplitude modes among the p-modes that we used in our phase 
analysis. Unfortunately, the highest amplitude periodicity is part of 
a triplet, that exhibits amplitude variation, and when the data are 
divided into subsets, a strong beating pattern in the phase analy- 
sis appears. There is a single periodicity fig with a relatively high 
amplitude. We have calculated amplitudes and phases of that peri- 
odicity in 30-day bins. The phases are stable to within 1 sec except 
for three points deviating by 5 a. The uncertainties of the phases are 
of the order of 0.5 sec. A similar flat trend was found using com- 
ponents of the quintuplet at 5287 ^iHz but with correspondingly 
larger uncertainties. It is beneficial to continue monitoring these 
periodicities for the remainder of the Kepler mission which should 
determine whether a companion (planet) exists. Spectroscopically, 
KIC 10139564 appears as a single star. 

5.5 Exotic species in the zoo 

Several regions in the amplitude spectrum appear to be frequency- 
crowded, unresolved, or "forked" which could indicate amplitude 
and/or phase variability. Investigating these regions has yielded 
unanticipated results, adding to the diversity of the pulsation zoo 
within KIC 10139564. 

The first exotic result is that of the regions around 5472 and 
5645 /jHz. While the region looks quite messy in the combined am- 
plitude spectrum (top panels of Figure|20l, the phases (shown in 
the bottom panels of Figurel20l. calculated using five-day sliding 
bins, drifted in a smooth fashion at a period of approximately one 
year. As such, we suspected an artifact related to the spacecraft but 
detected something even more interesting. The periodicities actu- 
ally occur above the Nyquist frequency so are reflected across it to 
the regions listed above. The true frequencies are centered around 
1 1521 and 1 1047 /iHz. Careful examination reveals these frequen- 
cies are, in fact, combinations from high-amplitude frequencies in 
the subNyquist region. For the periodicities near 11521 /iHz, there 
is a quintuplet formed from combinations of the main triplet (f33 
to fas) near 5760 /iHz. Similarly, the quintuplet near 1 1047 /iHz is 
caused by combinations of that same triplet (f33 to f3s) with the 
quintuplet near 5287 /iHz (fs to fi2). All combination frequencies 
we detected are listed in Table|3] While they are reflected across 
the Nyquist, they become smeared by the spacecraft's motion. The 
exposure cadence is fixed in the spacecraft reference frame. How- 
ever, the spacecraft is moving relative to the observed field, so the 
Nyquist frequency varies in the barycentric reference frame; this 
movement leads to a smearing of all frequencies whenever they are 
reflected across the Nyquist frequency. Therefore, these two messy 
regions are, in fact, reflections of combination frequencies above 
the Nyquist frequency. 

For other regions, this approach did not yield results, as shown 
in the right panels of Figure|20l Our best solutions produced an- 
other mixed set of exotic results. Amplitu de variability h as been 
previously observed in many sdBV stars ([Reed et al. 2 0071) . so it is 
not surprising also to encounter it here. One solution, therefore, is 
to investigate shorter sets of data over which the amplitudes and/or 
phases may be fairly stable. Our first approach was to divide the 
data into quarters, which is how the data are released to us. Fig- 
ure|2T]shows a region where this method resulted in a solution. The 
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Figure 20. The frequency regions near 5471, 5945 and 6003 /^Hz. The top panel shows the subNyquist regions, the middle panels the corresponding su- 
perNyquist regions and the bottom panel the phases of the subNyquist frequencies with errorbars. 



Table 3. The list of combination frequencies we found in the superNyquist 
region. We fitted these frequencies independently from those listed in Ta- 
ble|2] 



ID 


Frequency [/iHz] 


Period[s] 


Amplitude [ppt] 


S/N 


fci 


11046.7169(25) 


90.524634(21) 


0.042(8) 


5.1 


fC2 


11047.13862(42) 


90.5211779(34) 


0.251(8) 


30.2 


fC3 


11047.55681(34) 


90.5177513(28) 


0.310(8) 


37.2 


t"c4 


11048.3998(11) 


90.510845(9) 


0.093(8) 


90.6 


f C5 


11047.97594(14) 


90.5143173(11) 


0.754(8) 


11.2 


fC6 


11520.3240(6) 


86.8031145(45) 


0.175(8) 


21.1 


f C7 


11520.7520(7) 


86.799890(5) 


0.153(8) 


18.4 


tC8 


11521.1716(6) 


86.7967284(43) 


0.185(8) 


22.3 


fC9 


11521.5978(10) 


86.793517(7) 


0.110(8) 


13.2 


fcio 


11522.0152(8) 


86.790374(6) 


0.126(8) 


15.1 



region around 5840 /iHz was known to have portions of multiplets 
from previous quarters of data. Using Q5 through Q8 data, the am- 
plitude spectrum contained a complete septuplet above the detec- 
tion threshold. The peaks were not symmetric, indicating phase or 
amplitude variation, and when Q9 was added, several peaks dimin- 
ished below the detection limit. This sort of behavior is common 
for pulsations which are damped and re-excited. In Figure|2T] the 
septuplet frequencies of Q5 through Q8 are marked with dashed 
lines. The quarterly amplitude spectra reveal that no frequencies 
are above the detection limit during all five quarters, yet this region 



is quite frequency-rich. Quarters 5, 7 and 9 all show several sub- 
stantive peaks, which fit the previously established multiplet split- 
tings, and dotted lines in Figure[2T] indicate two previously unre- 
vealed frequencies. The average spacing between all these peaks is 
0.37 /iHz. The splittings are not identical, and there is likely some 
frequency wandering between quarters. This may be caused by am- 
plitude variations and/or phase interactions with other nearby peri- 
odicities which are also amplitude variable. We cannot claim all 
these periodicities are intrinsic to the star, but the evidence sug- 
gests this region contains an I mode. We also examined the 
multiplets near 3540 and 5049 /iHz by quarter, but no peaks were 
detected except those listed in Table[2] 

If quarterly data did not provide a solution for a given region, 
we examined them using shorter data sets. We calculated ampli- 
tude spectra for sets of data spanning 50, 60 and 70 days, avoid- 
ing the gap between days 270 and 290. The top panels of Fig- 
ure|22]and Table|4]show our solution for the region near 6001 /jHz. 
There appears to be some frequency wandering, but from these 
subsets of data, we estimate frequencies lie near 6001.10, 6001.47 
and 6001.88 /iHz, which give an overall separation of 0.39 /jHz, 
roughly the established rotational multiplet splitting. Most likely, 
this region contains a triplet with a separation of the spin period. 

The situation near 6173 [iHz is more complex (Figurel22l. 
largely because the amplitudes become very low between days 130 
to 220 and after day 330. Additionally, there are more frequencies 
involved in this region which covers ~ 3/aHz. It appears there is 
a triplet near 6170.83, 6171.17 and 6171.66 /jHz, which switches 
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Figure 21. The frequency region near 5840 fiHz. The top panel shows the 
entire 15-month data set, while subsequent panels show Q5 through Q9 
individually. The vertical dashed lines indicate frequencies detected using 
Q5 through Q8 data, and the dotted lines indicate two additional frequencies 
above the detection threshold in individual quarters that match the splitting. 
The horizontal lines denote the detection thresholds. 



power from the right pair to the center frequency to the outside 
pair prior to day 150. Then another triplet, which is modest at the 
beginning, gets even weaker in the middle but comes back at the 
end of the run. These frequencies are near 6172.39, 6172.83 and 
6173.32 fiHz, which would be another triplet. However, these can- 
not both be I = 1 triplets, since they are too close together, so one 
of them is likely a part of an I = 2 multiplet. This region, which is 
shown in the bottom panels of Figurel22l is very complex, and we 
can only roughly estimate the frequencies of Table|4] Any interpre- 
tations are fairly speculative. We estimate that this region contains 
two distinct triplets with splittings of the rotation period, but we 
certainly cannot rule out a single multiplet. 

The last region we discuss is near 6235 /jHz (Figurel23V In 
this region, there is a high amplitude peak in the overall amplitude 
spectrum with two sidelobes, none of which prewhiten well. Break- 
ing the data into 60-day subsets to resolve the frequencies clarifies 
what is occurring. In this case the main frequency f42 at 6234.7 is 
frequency and phase stable with an amplitude that diminishes after 
day ~80. There are two other periodicities, both with frequencies 
and amplitudes that vary during the observations. The higher fre- 
quency (f43) drifts from 6235.07 to 6235.02 /jHz and back again 
during the run, while the lower frequency drifts from 6234.4 to 
6234.2 /jHz prior to day 180, after which it is no longer detectable. 
This is shown in Fi gure|23l where the original amplitude spectra are 
shown in the top panels, the amplitude spectra prewhitened by £42 



Table 4. Frequencies detected from analyses discussed in i|5.5l This list also 
includes prewhitened frequencies from Table|2] 



ID 


Frequency [^Hz] 


Au 


I 


m 


Xi 


5838.64 


X2-0.33 




-4 


x 2 


5838.97 


f 3 6-0.41 




-3 


f36 


5839.380 


X3-0.31 




-2 


X3 


5839.694 


X4-0.41 




-1 


A4 


JQ4U.1UU 




A 

4 


U 


x 5 


5840.479 


X 4 +0.38 




+ 1 


x 6 


5840.841 


Xs+0.36 




+2 


Xr 


5841.198 


X 6 +0.36 




+3 


f 3 7 


5841.577 


X7+O.38 




+4 


fiS 


6001.10 


f 3 9-0.37 




-1 


f39 


6001.47 


- 


1 





x 8 


6001.88 


f 3 9+0.41 




+ 1 


x 9 


6057.4 


- 


- 




X10 


6057.6 


- 


- 




Xn 


6170.83 


X12-O.34 




-1 


X12 


6171.17 




1 





X13 


6171.66 


X12+O.49 




+1 


X14 


6172.39 


Xis-0.44 




-1 


X15 


6172.83 




1 





X16 


6173.32 


Xi5+049 




+1 


X17 


6234.3 


f 47 -0.39 




-1 


f47 


6234.7 




1 





f48 


6235.05 


f 47 +0.32 




+1 



in the middle panels and the amplitudes and phases in the bottom 
panels. Clearly the amplitudes of all frequencies are variable. The 
amplitudes for f 4 g are only significant for the fits centered at days 
90 and 150 while the amplitude for f 4 7 continually declines until 
by day 380 it is no longer significant. 

We also attempted to disentangle the region near 6314 /jHz us- 
ing these techniques, but we were not able to establish reasonably 
stable frequencies. The estimated frequencies for these regions ap- 
pear in Table|4] 
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Figure 22. Amplitude spectra for the regions around 6001 and 6173 /iHz. Left panels show the entire data set, while the remaining panels show 50-day 
segments as labelled. Dotted lines indicate likely frequencies listed in Table|4] 



6 SUMMARY 

KIC 10139564 is the only p-mode dominated sdBV star observed 
by Kepler. It was discovered during the survey phase of the mis- 
sion and has been observed continuously in SC mode since the 
start of the second year. The SC mode allows us to study short 
period pulsations that would be unobservable with LC sampling. 
KIC 101539564 contains frequencies higher than the Nyquist fre- 
quency associated with the SC sampling interval, but thanks to the 
long duration of these observations the barycentric correction in- 
duces frequency modulations that resolves the ambiguity of modes 
reflected across the Nyquist frequency into the subNyquist region. 

The amplitude spectrum calculated from 15 months of nearly 
continuous data shows more than 50 periodicities in a pulsation 
zoo, which includes both anticipated and exotic species. We de- 
tected the anticipated p-modes, found several g-modes and two 
periodicities between the p- and g-mode regions. These pos- 
sibly mixed modes have been observed in some other sdBV 
stars dReed et all l20ld : iBaran et al.l 1201 lah . and are particularly 
abunda nt in the complex pulsator 2M 193 8+4603 d0stensen et al.l 
l2010ch . The highest amplitude is around 8 ppt, and only a few other 
frequencies are of comparable height. Most of the signal is below 
0.5 ppt, which would be difficult to detect from the ground. 

Almost all periodicities are split into frequency multiplets. We 
interpret this splitting as caused by rotation and derive an average 
rotation period of 25.6±1.8 days. The observed splittings vary by a 
small amount, leading to an uncertainty in the rotation period. This 
is most likely caused by small differences i n the Ledoux constant a s 
has been observed in models calculated bv lCharpinet et alj J2002h . 
The scale of variation in the Ledoux constant matches the variation 
in frequency spacings. The exact rotation rate can be established 
once the order n and degree I are identified to better constrain the 
Ledoux constant. Most gravity modes appear as triplets, and their 
frequency differences are consistent with identifying them as I = 1 
modes. 



They also show nearly equal period spacings, which we inter- 
pret as asymptotic overtones of a single degree. Multiplets indicate 
they are I = 1 modes, but the period spacings are too large to be of 
consecutive overtones. They are more likely 207 or 310 s which is 
dissimilar to the 250 s observed fo r the other Kepler-observed pul- 
sating sdB stars dReed et al.ll201 1). However, those other stars are 
all g-mode dominated and so are cooler than KIC 10139564. The 
only other sdB star that is a p-mode dominated hybrid is Bal09. 
Ground-based observations detected 20 g-mode periodicities yet 
there is no indication of consistent period spacings between them. 
So at this time Bal09 is the only pulsating sdB star with rich g- 
mode spectrum which does not show asymptotic g-mode period 
spacings. It is also both cooler and more evolved (lower log g) than 
KIC 10139564 and so may not be any more suitable for comparison 
than the other g-mode pulsators. Theoretical works in progress (by 
VVG) are seeking to explain this discrepancy. 

Several of the multiplets show signal excesses at the frequen- 
cies of likely missed components and are indicated in the figures 
with vertical dash-dotted lines. We included them in our list as the 
S/N = 4 threshold is adopted to separate intrinsic peaks from ran- 
dom noise. When there is good reason to anticipate pulsation fre- 
quencies, in this case multiplet splitting, the significance threshold 
does not necessarily apply. We do not claim that any of these fre- 
quencies are intrinsic to the star, yet we include this speculation to 
compel future investigation with additional data. 

Recent spectroscopic measurements indicate KIC 10139564 is 
not a radial velocity variable and allowed us to refine its tempera- 
ture and gravity, placing it well away from t he typical hybrid pul- 
sator region (see Fig. 2 of 0stensen 2010b). This raises the pos- 
sibility that low-amplitude, g-mode pulsations are common place 
in V361 Hya class stars, but an instrument with Kepler's accuracy 
is required to observe them. A re-examination of pulsation driv- 
ing models woul d seem necessary as the most recent g-mode driv- 
ing publications jjefferv & Saiol2006tlHu et alj201 1) cannot drive 
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Figure 23. The frequency region near 6235 fiHz. Original amplitude spectra (top), residuals after prewhitening the strongest frequency (middle) and the 
amplitudes (bottom left) and phases (bottom right) for the frequencies at 6234.7 (filled circles) and 6235.0 fiHz (open squares). 



pulsations to the temperature of KIC 10139564. Improvements for 
diffusion of iron-group elements (in progress by HH) seems likely 
to solve this issue. 

Among the p-modes at high frequencies, we have noticed sev- 
eral regions to be unresolved in frequencies, and our first guess led 
us towards stochastic oscillations as a possible explanation. De- 
tailed analysis showed, however, that they are likely not stochastic 
in nature, but rather require multiple explanations. The regions all 
appear as multiplets and frequency splittings are consistent with 
our derived rotation period. Two regions are caused by reflections 
across the Nyquist frequency of combination frequencies, which 
are smeared by the difference in spacecraft and barycentric refer- 
ence frames. Others are caused by highly variable amplitudes, and 
some appear to be changing their frequencies slightly. In the com- 
bined data set, this leads to multiple-peaked or even transient peri- 
odicities with amplitudes that are greatly reduced. 

As a result of the abundant multiplets, we were able to as- 
sociate most frequencies with modes, which is a very rare result 
for sdB V stars. Also of interest are the probable detections of I = 3 
or greater modes. In many other sdBV stars, freque ncy density 
has st rongly suggested that higher I modes are excited I IReed et alj 
1 20071) , yet because of geometric cancellation, there has been some 
debate whether we would most likely observe 2 = 3 or 4 modes. 
Mode density in t his star also creates a new problem because cur- 
rent models (e.g. Charpinet et al. 20o3) cannot produce as many 
modes as are observed. From a sample of representative models, the 
overtone spacings for I = to 4 averaged from 904 to 939 /iHz. Near 
the fundamental order, this could be reduced to about 700 fiHz. For 
the better-established sequences, we detect overtone separations of 
237, 124, 176, 241 and 234/iHz. Such a dense spectrum is clearly 
challenging for current models yet could go a long way to explain 
the high-density observed in several sdBV stars. KIC 10139564 
seems well placed to solve these issues, and the many multiplets 



(especially those of higher degree) should provide constraints on 
the rotation axis as well. 

Only one high-amplitude p-mode is detected to be a single 
peak. Its amplitude is almost 4ppt. We used this peak to match 
amplitudes between months to reduce residual signal as much as 
possible. This helped to diminish residuals for most of the peaks, 
which could mean either the change of intrinsic amplitudes of those 
periodicities is consistent or they are stable and any difference in 
amplitude between quarters is caused by systematics. 

Fifteen months of continuous data coverage opened a way to 
study phase stability of detected periodicities to search for possible 
companions to KIC 101039564. Although we have detected more 
than 50 periodicities, we could not use most of them for checking 
phase stability. Close neighbouring peaks or residuals, which in- 
dicate amplitude/phase variation, cause a beating pattern in both 
amplitude and phase. We found only one well-suited frequency 
around 5472 /iHz, which we used to look for phase stability. It 
is still affected by some nearby low-amplitude, complex signals. 
However, the beating between the peak and that complex signal is 
about 10 days, so the beat period is shorter than the period from 
rotational splitting. As a result of our phase analysis, we have not 
noticed any periodic variation in phase that could be interpreted as a 
companion to KIC 10139564. There seems to be a real deviation in 
phase between day 420 and 520, but we are not yet able to inteipret 
this variation. 

The results we presented in this paper clearly show 
KIC 10139564 to be very interesting and allowed us to detect fea- 
tures not accessible from ground-based data. As most frequency 
regions contained multiplets, mode identifications could be deter- 
mined for most of the peaks we detected. The identified modes 
should pose strong constraints on stellar models and should, ob- 
viously, allow the exclusion of non-unique solutions, making the 
theoretical study of this object easier. In fact, this is the only sdBV 
star with so many modes identified. As a natural implication of the 
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pulsation variety and all the remarkable results we obtained from 
only 15 months of data, we definitely should continue monitoring 
this star to confirm several features as well as to uncover more prop- 
erties still hidden in the data. 
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